Theory of excitonic spectra and entanglement engineering in dot molecules 
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We present results of correlated pseudopotential calculations of an exciton in a pair of vertically stacked 
InGaAs/GaAs dots. Competing effects of strain, geometry, and band mixing lead to many unexpected 
features missing in contemporary models. The first four excitonic states are all optically active at small 
interdot separation, due to the broken symmetry of the single-particle states. We quantify the degree of 
entanglement of the exciton wavefunctions and show its sensitivity to interdot separation. We suggest 
ways to spectroscopically identify and maximize the entanglement of exciton states. 
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The small size of semiconductor quantum dots [1] drives 
speculations that they may provide a physical representa- 
tion of a quantum bit (qubit) that supports a superposition 
of "0" and "1" HUBBlSIl- Some proposed representa- 
tions of qubits include electron spin states [ 2] and the pres- 
ence or absence of an electron, hole, or electron-hole pair 
fexciton") . m. lol Hol nm Registers of qbits might be realized 
in coupled quantum dots, self-assembled by strain-driven 
islanding of InGaAs on a GaAs substrate 1 12]. In one pos- 
sibility, an electron represents qubit A and a hole represents 
qubit B, while the qubit states are the occupation of either 
the top (T) or bottom {B) dot Li,^ia,JJJ. This quantum 
register must store entangled states. Predicting entangle- 
ment requires a theory of the electronic structure of the dot- 
molecule, including single-particle and correlation effects. 
Most modeling of dot-molecules has been done in siiigle- 
band effective-mass approximation [7, 13, 14, T?, 16,T7|. 
For two equivalent dots, this treatment leads to single- 
particle electron and hole orbitals forming bonding and an- 
tibonding combinations: 






where ej (bb) represents an electron in the top (bottom) dot 
and h-Y (/ib) represents a hole in the top (bottom) dot. In 
this picture, as in an HJ molecule, the single-particle bond- 
ing state energy E^[d) decreases as the interdot distance d 
decreases while the bonding hole state energy E^{d) in- 
creases. Simple direct products of single-particle states, 
e.g. \(t)\<j)l) = i[|eT/iT> + lefi/iB) + Ict/ib) + leB/ix)], 
are unentangled. In contrast, the desirable maximally en- 
tangled Bell states superpose either exciton or dissociated 
states, but not both: 

\a) = -^{\e-Bh-Q) + [ct/it)); exciton, bonding, 
\h) = -7|(|eB/iT) + Icx/iB)); dissociated, bonding, 
\c) = -7|(|eB/iT) — Ict^b)); dissociated, antibonding, 
\d) = -^{\e^h^) — \ejhj)); exciton, antibonding. (2) 



Recent optical experiments on vertically-stacked double 
dots 10, IiHi claimed to show entangled excitonic states, 
but the evidence for entanglement is indirect and based 
on a symmetric model underlying Eq. Q. Unlike simple 
symmetric molecules like HJ, double-dot "molecules" of 
stacked InGaAs dots are made of -^10^ atoms and have 
complicated interactions such as alloy fluctuations, strain, 
multi-band (e.g. light-heavy hole), inter-valley (T-X), and 
spin-orbit couplings not included in the symmetric molec- 
ular case. To properly simulate these double dots, we have 
performed detailed atomistic pseudopotential calculations, 
including correlation, on a realistic dot molecule. In this 
Letter we report on new insights into the exciton state: all 
states are optically active at short distances, entanglement 
is small except at a critical dot separation dc at which the 
low energy exciton is darkened, yielding a spectroscopic 
signature of entanglement. 

To understand previous theoretical treatments of exci- 
tons in a dot molecule and to set a reference to which our 
atomistic results will be compared, we describe a generic 
two-site tight-binding Hamiltonian in the basis of products 
of electron and hole single-particle states Ict/it); leB^iT); 
lex/iB), Icb/ib): 
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Here {ej,eg ,ej[,ej^} are the electron and hole on-site 
energies, {tejih} are the tunneling matrix elements, and 
{W^ ,U]^ ,Uf^ ,Uf^} are the electron-hole Coulomb 
matrix elements. The extra parameter IS.E that will be 
used later in entropy discussion; initiall y w e set Ai? = 0. 
A simplification, followed in Refs. H Efl, Hil l, is to set: 



ie = ^h = t, ^J = £? and el = e^, with intra-dot 



TT 



Coulomb energies U^^ 



ttBB 
^eh 



-h' 



U, and neglecting 



inter-dot terms C/™ and C/^^. With this simplification, the 




4 6 8 10 12 14 16 18 20 22 
base-to-base separation (nm) 

FIG. 1 : Exciton energies as a function of the interdot separation 
for (a) our pseudopotential many-body results, (b) a model cal- 
culation A_B = ^/Jlij^TU^. 



Hamiltonian of Eq. (|3j yields, in increasing order of en- 
ergy, the four e-h states: 



|l)=(|a)-7i|&))/VrT^, 
|2) =|c), |3) = \d), 



l4)=(|a)-72|&))/7rT^, 
7i,2 = \U± V(4t)2 + UA /4t 



(4) 



We see that in the simplified model states |2) and |3) are 
fully entangled pure Bell states [viz. Eq. (jjjl] that are spa- 
tially antisymmetric (anti -bonding) and therefore optically 
dark. In contrast, states |1) and |4) are not fully entan- 
gled and have some symmetric (bonding) character mak- 
ing them optically allowed. Assuming that the tunneling 
integral t{d) decays with interdot separation d, the simple 
model gives the level order shown in Fig.^^). In the simple 
model, the exciton |1) shifts to the red as d decreases, and 
the separation of the two bright states |1) and |4) increases 
as Aii^ = y^(4t)2 + If^. Furthermore, the order of the 
levelsis |l)=bright^ |2)=dark^ |3)=dark^ |4)=bright. 
This is in apparent agreement with exp eriments that show 
the same qualitative behavior ISL flOLflll l. spurring hope that 
the theoretically predicted high degree of entanglement in 
this system could be experimentally realized to the benefit 
of quantum computing. 

There are reasons to doubt the simple homonuclear 
diatomic-like analogue of dot molecules of Eq. Q and 
Fig. [Hb). Actual self-assembled quantum dots contain 
~10^ atoms, while the dots themselves are strained by 
the host matrix and subjected to random alloy fluctuations. 
Thus, a "molecule" made of two dots does not necessarily 
behaves like homonuclear H2, but could behave like a het- 
eronuclear molecule (e.g. HE) since strain and alloy fluc- 
tuations distinguish the dots, e^ ^ e^. Furthermore, the 
electronic properties of such dots cannot 1 18] be accurately 
modeled by simple single -band effective-mass approaches: 



CO 



CD 



o 



1 1 


(b) 


GaAs 


d = 22.6 nm 


1 1 


ms^'"^ 


d = 17.0 nm 


d 


,10 nm^ 

/R8i^|2nm 


II 


d = 11.3nm 


II 




' 12nm ' 


d = 10.2 nm 


1 


d = 9.8 nm 




; 1x2 


d = dc = 8.5 nm 


i ;x2 


d = 7.9 nm 


.1 


d = 7.3 nm 


d = 6.8 nm 


d = 5.7 nm 
, 1 


d = 5.1 nm 
■ 1 



1.24 1.26 1.28 

transition energy (eV) 



1.30 



FIG. 2: (a) Emission spectra in a pair of vertically stacked In- 
GaAs/GaAs dots, (b) Dot geometry, including a two monolayer 
(0.56 nm) InGaAs wetting layer and graded composition profile. 



coupling between a large number of bands alters electron 
and hole localization, changing the Coulomb matrix ele- 
ments. Finally, the assumption of equal tunneling for elec- 
tron and hole, tg = th, is questionable given the large mass 
ratio, rrie/rnhh ~ 1/6, of electrons and heavy holes in the 
GaAs barrier between the dots. In fact, we see below that 
band mixing even changes the sign of th- Thus, a more 
complete theoretical treatment is warranted. 

We simulate the InGaAs/GaAs dot molecule at a range 
of inter-dot spacings, using a computational approach 
that successfully describes single InGaAs/GaAs dots [19]. 
Specifically, we describe the single-particle properties with 
an atomistic empirical pseudopotential method, with the 
wavefunctions expanded in a set of Bloch states of the con- 
stituent materials over many bands and wave vectors 1201 . 
The theory includes multi-band and multi-valley coupling, 
spin-orbit interaction, and anisotropic strain effects. To cal- 
culate correlated e-h states, we include excitonic effects 
in a low-order configuration interaction expansion, as in 
Ref. I21I, calculating all Coulomb and exchange integrals 
explicitly from the single-particle wavefunctions. The dot 
geometry has been chosen to resemble the experimental 
system studied by Bayer et al [8]. As shown in Fig.|3b), 
the dots are 12 nm x 2 nm truncated-cone-shaped, with a 
linear composition gradient varying from Ino.5Gao.5As at 
their bases to pure InAs at their tops. Fig. |3] shows the 
calculated single -particle energies and wavefunctions. We 
plot our calculated correlated e-h energies. Fig. ^a) and 
corresponding absorption spectra, Fig.|2l 

By projecting our numerically calculated correlated e-h 
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FIG. 3: Single-particle energies as a function of interdot separa- 
tion and single-particle state densities from our pseudopotential 
calculations. The two translucent isosurfaces enclose 75% and 
40% of the total state densities. The physical dot dimension are 
shown in grey, with base-to-base separation d=5.1 nm. 



energies vs. d onto a generalized form of H, Eq. (O, we 
have determined the effective distance-dependent Hamil- 
tonian parameters, shown in Fig. |3 These represent re- 
alistic values for the simplified model parameters contem- 
plated in Ref. 2. Inspection of the parameters from Fig.|3 
and the comparison within Fig.^show that our results dif- 
fer in crucial ways from the simple assumptions made in 
Refs. ^, U^ [ni We next discuss the salient features of 
these differences and their physical implications. 

(i) The energies of exciton |1) and |2) blue-shift as d 
decreases, in contrast with the red-shift expected from the 
simple model. Two effects are responsible for this sur- 
prising blue-shift, (a) From a single-particle view. Fig. |3l 
shows that as the interdot separation decreases, the en- 
ergy of both hole states /iq and hi move to lower val- 
ues, while the molecular bonding-antibonding picture of 
Eq. O would predict that the bonding /iq level will moves 
to higher energies. The downward shift of the single- 
particle hole level /iq with decreasing d contributes to the 
upward shift of the lowest excitons |1) and |2) observed 
in Fig. [Ha). The reason that the single-particle states /iq, 
hi move to lower energies as d decreases is their symmetry 
broken character: /iq is localized on B and hi on T, as seen 
in the densities, Fig.|3] This localization is reflected in the 
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FIG. 4: Effective parameters for the two-site Hamiltonian, 
Eq. J3j, distilled from our pseudopotential calculations. 



small tunneling matrix element for holes th, in Fig-S and is 
due to the heavier hole mass and the strong strain-induced 
potential barrier |22] between the dots. In contrast, the 
light-mass electrons have a large tunneling matrix element 
te and follow the bonding-antibonding picture of Eq. Q, as 
evidenced by the calculated density. Fig. |3l exhibiting de- 
localization on both dots, (b) From an interacting-particle 
view, the blue shift of |1) and |2) with decreasing d is 
caused by the decrease in the Coulomb elements U^ and 
U^^ with reduced interdot separation (shown in Fig.Uc)), 
due to delocalization of the exciton on both dots. 

(ii) At large d > \{) r\m the order of excitons |1), |2), |3) 
and 1 4) is bright, bright, dark and dark, in contrast with 
the simple model predicting the order bright, dark, dark, 
bright. The large d behavior of our pseudopotential cal- 
culations can be understood in the tight-binding language: 



differences in on-site energies (ej — ej^) — (e^ — £h) are 
greater than hopping elements tg and th> (Fig- 0- With 
these assumptions, the exciton states, in increasing order 
of energy, are given by |1) = Ict/it): |2) = \e^h^), 
|3) = Ict/ib), |4) = lefi/iT)- States |1) and |2) are bright 
since they are symmetric and have large e-h overlap. In 
contrast, states |3) and |4) are dark since they are not sym- 
metric and have low e-h overlap. These four eigenstates 
are obviously not entangled, while the simple model pre- 
dicts full entanglement. 

(Hi) Exciton [1) and |2) anticross at the critical distance 
dc at which point |1) becomes dark. However, all exci- 
tons are bright at d < dc, in contradiction with the simple 
model. This can be understood as follows. At dc (—8.5 
nm for our specific case) the basis states \ejhi) and {e^h^) 
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FIG. 5: Degree of entanglement after Eq. J5}- The dashed Une 
indicates the value AE = el — £?=2.9 meV from our simulations. 



nm arises from our intrinsic dot energy separation Ai? = 
el - £^=2.9 meV. 

Extracting the entanglement from this exciton state for 
use in quantum computing may require the separation of 
the electron and hole while maintaining phase coherency. 
This experimental challenge might be accomplished by 
driving the particles to nearby dot molecules using an in- 
plane electric field. 

In conclusion, we find that the entanglement entropy 
reaches a maximum value (of 80% in our case) at a critical 
interdot separation and decays abruptly to zero at smaller 
and larger separations. We suggest that the distance d,. can 
be identified using photoluminescence spectra, by noting 
two closely spaced exciton peaks with a darker lower en- 
ergy peak. 

Work supported by the DOE SC-BES-DMS Grant No. 
DE-AC36-99GO10337. 



are nearly degenerate. Now hopping elements t^, and th 
will split this near degeneracy into symmetric and antisym- 
metric combinations |2) = \a) and |1) = |6), respectively. 
Whether the ground state is symmetric or antisymmetric is 
decided by the respective signs of tf. and th- If they have 
opposite signs, the symmetric e-h state has lower energy. 
This is the case in the single band effective mass approx- 
imation lll6ll where the single-particle hole and electron 
states have pure S-envelope-character creating Vssa bonds 
12311 with positive th and negative t^. However, in a re- 
alistic case, the presence of inter-band mixing 1 18] (S-P as 
well as heavy and light hole) leads to a V^p^-like hole bond, 
with th < 0. In fact, the confining potential for holes at- 
tracts the light-hole P-like component of the hole states to 
the bonding region in-between the dots. If t^ and t^ have 
the same sign (viz, our Fig. |4t)), the antisymmetric state 
(dark) is below the symmetric (bright) state . 

(iv) The degree of entanglement reaches a maximum at 
a critical interdot separation dc-' Since |1) and |2) ap- 
proach Bell-states \a) and \d) at d = dc, we expect high 
entanglement. We quantify the degree of entanglement di- 
rectly from our pseudopotential e-h density matrix, using 
the von Neumann entropy of entanglement |24], 



5 = - Tr pe loga Pe 



TrphlogaPh, 



(5) 



where pe and ph are the reduced density matrices of the 
electron and hole, respectively. For exciton |1) we find 
S = 0% for d > 10 nm and a pronounced peak at dc 
with S = 80%. dc is determined by a balance between 
interdot strain coupling and intrinsic dot energy differences 
(alloy fluctuations, in our calculations). To generalize our 
results to a class of dots we allow in Eq. (|3j a generic 
fluctuation AE. We have calculated S{d, AE), as shown 
in Fig. |5] using our fitted model Hamiltonian, Eq. (|3j and 
Fig. |4] This shows how various degrees of entanglement 
can be engineered. We note that the specific case dc = 8.5 



[1] D. Bimberg, M. Grundman, and N. N. Ledentsov, Quan- 
tum Dots Heterostructures (John Wiley & Sons, New York, 
1999). 
[2] D. Loss and D. P. DiVincenzo, Phys. Rev. A 57, 120 (1998). 
[3] P Zanardi and F. Rossi, Phys. Rev. Lett. 81, 4752 (1998). 
[4] A. Imamoglu et al. Phys. Rev. Lett. 83, 4204 (1999). 
[5] P Zanardi and F Rossi, Phys. Rev. B 59, 8170 (1999). 
[6] E. Biolatti, R. C. lotti, P. Zanardi, and F Rossi, Phys. Rev. 

Lett. 85, 5647 (2000). 
[7] F Troiani, E. Molinari, and U. Hohenester, Phys. Rev. Lett. 

90, 206802 (2003). 
[8] M. Bayer et al. Science 291, 451 (2001). 
[9] M. Korkusinski and P Hawrylak, Phys. Rev. B 63, 195311 
(2001). 
[10] M. Korkusinski et al. Physica E 13, 610 (2002). 
[11] M. Bayer Physica E 12, 900 (2002). 
[12] I. Shtrichman et al, Phys. Rev. B 65, 081303 (2002). 
[13] L. R. C. Fonseca, J. L Jimenez and J. P. Leburton Phys. Rev. 

B 58,9955 (1998); 
[14] C. Yannouleas and U. Landman, Phys. Rev. Lett. 82, 5325 

(1999; ibid. 85, 2220 (2000). 
[15] M. Koskinen, S. M. Reimann, and M. Manninen, Phys. Rev. 

Lett. 90, 066802 (2003). 
[16] F Troiani, U. Hohenester, and E. Molinari, Phys. Rev. B 65, 

161301 (2002). 
[17] T. Fujisawa, D. G. Austing, Y. Tokura, Y Hirayama, and 

S. Tarucha, Phys. Rev. Lett. 88, 236802 (2002). 
[18] L.-W. Wang et al, Appl. Phys. Lett. 76, 339 (2000). 
[19] A. J. Williamson, A. Franceschetti, and A. Zunger, Euro- 
phys. Lett. 53, 59 (2001); J. Shumway et al, Phys. Rev. B 
64, 125302 (2001); G. Bester, S. Nair, and A. Zunger, Phys. 
Rev. B 67, 161306 (2003). 
[20] L. W. Wang and A. Zunger, Phys. Rev. B 51, 17398 (1995); 

Phys. Rev. B 59, 15806 (1999). 
[21] A. Franceschetti et al, Phys. Rev. B 60, 1819 (1999). 
[22] W. D. Sheng and J. P Leburton, Appl. Phys. Lett. 81, 4449 

(2002). 
[23] W. A. Harrison, Electronic Structure and the Properties of 



Solids (Freeman San Francisco, San Francisco, 1980). Phys. Rev. A 53, 2046 (1996). 

[24] C. Bennett, H. Bernstein, S. Popescu, and B. Shumacher, 



